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Abstract 



We propose the construction of highly directional emission microlasers using two-dimensional high-index semiconductor waveg- 
uides as open resonators. The prototype waveguide is formed by two collinear leads connected to a cavity of certain shape. The 
proposed lasing mechanism requires that the shape of the cavity yield mixed chaotic ray dynamics so as to have the appropiate 
(phase space) resonance islands. These islands allow, via Heisenberg's uncertainty principle, the appearance of quasi bound 
states (QBS) which, in turn, propitiate the lasing mechanism. The energy values of the QBS are found through the solution of 
the Helmholtz equation. We use classical ray dynamics to predict the direction and intensity of the lasing produced by such 
open resonators for typical values of the index of refraction. 

PACS: 42.55.Sa, 05.45.-a 
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Studies on microcavities which serve as optical res- 
onators have led to the design of semiconductor micro- 
lasers, see for example [1]. One way to obtain high 
performance microcavities is through the excitation of 
Wispering- Gallery Modes (WGM), first observed in 1961 
using a spherical sample of CaF 2 : Sm 2+ [2]. WGM 
living in microcavities with spherical (or circular) geom- 
etry were studied in polyestyrene spheres [3], spherical 
doplets [4] , spheres of fused quartz [5] , and in semiconduc- 
tor disks [6] . In Ref. [2] it was also suggested that WGM 
could be tuned by deforming the cavity from a spherical 
shape. This idea has been explored in spheroidal doplets 
[7] , fused silica spheroids [8] , deformed fused silica spheres 
[9,10], deformed disks [11], and other geometries [12]. The 
WGM in a deformed disc were first analyzed in terms of 
chaos theory in Ref. [13], and later a similar analysis was 
made for deformed doplets [14]. Furthermore, based on 
features of nonlinear dynamics, the asymetric resonant 
cavity (ARC) as a resonator for semiconductor lasers was 
proposed [15]. The ARCs are two-dimensional (2D) cir- 
cular resonators with quadrupolar deformations. Small 
deformations from circular geometry were shown to spoil 
the Q of WGM but it was found that for higher defor- 
mations &ow-iie-shaped resonances, instead of WGM, are 
responsible for the improvement in the power and direc- 
tionality of the laser cmmission [16]. 

Here we propose the use of a 2D waveguide system as 
an open resonator with similar topological properties as 
those of the ACRs with large deformations. We believe 
that the great advantage of an open resonator may be 
in that no pedestals or couplers (pumpers) close to the 
cavity are needed, as is the case with the currently inves- 
tigated 3D and 2D microlasers, see e.g., Refs. [8,6,11,16]. 

The waveguide system is formed by a cavity connected 
to two collinear semi-infinite leads of width d extended 



along the x-axis. Our prototype waveguide has the ge- 
ometry of the so-called cosine billiard extensively studied 
in [17-21]: it has a flat wall at y — and a deformed wall 
given by y(x) = d + a[l — cos(2ttx/L)], where a is the 
amplitude of the deformation, and L is the length of the 
cavity. 

We shall now discuss dynamical features of this system 
using, firstly, ray optics since it is assumed that the wave- 
length of light A is much smaller than the amplitude a of 
the cosine profile. 

The set of parameters (d,a,L) = (1.0,0.305,5.55) 
is known to produce mixed chaotic dynamics [19-21] 
(an incomplete ternary horseshoe [19]). The complete 
panorama of the ray dynamics is given by the Poincare 
Map (PM), which is the plot of the orbits of a set of rep- 
resentative initial conditions as these cross, in a specified 
direction, a given surface of section in phase space [22]. 
For our purposes the surface of section can be either the 
botton (flat) or the top (cosine) boundary. Thus, each 
time a ray impinges on the top wall we plot the position x 
and the sine of x, where \ is the angle the ray makes with 
the normal to the top boundary, see Fig. 1(a). Similarly 
we obtain another plot, see Fig. 1(b), when we use the 
bottom wall as a surface of section. Small variations in 
the geometrical parameters yield qualitatively the same 
topology in phase space. In Fig. 1(a) we clearly identify 
three islands: one large at sin(x) = and two smaller 
ones centered at sin(x) = 0.3. Motion within these is- 
lands correspond to bounded motion inside the cavity 
(i.e. not accesible to scattering trajectories) and are pro- 
duced by rays bouncing in the neighborhood of stable 
period-one and period-four fixed points, respectively [23]. 
With the surface of section at the top boundary, the phase 
space topology of this cavity is very similar to that of the 
ARC with large deformations [16]. Observe that even 



though both figures, 1(a) and 1(b), are constructed us- 
ing the same orbits in position space, their phase space 
topology are different. As shown below, we shall take 
advantage of this difference. 
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FIG. 1 Poincare Maps for the closed cavity. The Poincare 
surface of section is located at (a) the deformed wall: 
y(x) — d + a[l — cos(2ttx/L)], and (b) the flat wall: y(x) = 0. 
The parameters (d,a,L) = (1.0,0.305,5.55) were used. Some 
critical lines for refractive scape are also shown. 

To study the wave scattering phenomena we numeri- 
cally solve the Schrodinger equation for our waveguide, 
as in [18,19], since we know [24,25] that the problem of a 
TM wave inside a 2D waveguide with Dirichlet boundary 
conditions (Hclmholtz equation) is equivalent to that of 
a quantum wave in a 2D billiard (Schrodinger equation). 

It has been shown [19,20] for this waveguide with ge- 
ometrical parameters (d,a,L) = (1.0,0.305,5.55) that a 
plot of the Landauer conductance G presents sharp reso- 
nances at energies where the wave functions either, tun- 
nel into classical resonance islands (forbiden phase space 
regions for classical scattering) , or have their support ar- 
round them. Those scattering wave functions that tunnel 
into classical islands are very similar to a set of eigen- 
functions of the corresponding closed cavity [21], hence 
we shall refer to them as quasi bound states (QBS). It 
was also found [21] that QBS may be of two types only: 
M-type or I-type, see Fig. 2. The M-type (I-type) corre- 
sponds to a wave function with support on the period-four 
(period-one) resonance. These patterns are produced by 
waves incident from the left, modes 3 and 6 respectively. 
Note that there is no evidence of reflection, which is a 
general feature of the lowest modes, as shown in Refs. 
[18,19]. 

In Fig. 3 we present a plot of G for a range of energies 



supporting 20 propagating modes. The insets (a) and (b) 
enlarge two sharp dips of G corresponding to the QBS of 
Fig. 2(a) and 2(b), respectively. 




FIG. 2 Wave functions | ^(x,y) | 2 calculated at resonance, 
(a) E = 4008.8583, (b) E = 4027.0597. These QBS corre- 
spond to modes (a) 3, and (b) 6. The cavity is defined in the 
range x = [0,5.55]. E = El/E* ~ M^ ax Tv 2 ; where El is the 
expression for the energy in the leads, E* = ft 2 /(2m e d 2 ), and 
Mmax is the largest transversal mode m beyond which the 
longitudinal wave vector k m = \J (2m e E L /h 2 ) — (m 2 -K 2 /d 2 ) 
becomes complex. 




FIG. 3 Dimensionless conductance G as a function of E. 
Insets (a) and (b) enlarge the dips whose energies where used 
to calculate the wave functions of Figs. 2(a) and 2(b), respec- 
tively. 

If the waveguide system is constructed with a semi- 
conductor material having a refraction index n and we 
introduce a wave into one of the leads with an energy 
corresponding to a sharp resonance in G, the cavity acts 
as an open resonator since the light will bounce in the cav- 
ity region following an M-type or an I-type orbit. Also, 
we assume silvered surfaces along the leads so that light 
inside the waveguide undergoes especular reflections at 
the silvered surfaces and have a chance to escape only 
in the cavity region. According to Snell's law there is a 
critical angle, \c = sin _1 (l/n) such that, if \ > Xc the 
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ray light reflects specularly, while if \ < Xc the ray light 
refracts. In Fig. 1 some critical lines for refractive escape 
sin(xc) are ploted for different refractive indices. For a 
given value of n, orbits corresponding to points above 
(below) the critical line will reflect (escape). 

Since the QBS form M-type or I-type patterns we can 
use ray dynamics to estimate the direction (using Snell's 
law) and the intensity (using Fresnel transmitance and 
rcfractancc) of the emission for a given refractive index 
n. We do this by following the motion of an ensamble of 
rays with initial conditions inside the resonance islands 
and by counting the number of times the rays visit each 
of the elements of a mesh defined both, inside and out- 
side (to consider escape from) the cavity. The mesh is 
a 555 x 555 square lattice which spans the whole plots 
of Figure 4. At the end of the counting each mesh ele- 
ment with coordinates (x, y) has a weight p(x, y) which is 
the classical counterpart of the wave function | &(x, y) | 
[26]. In Fig. 4 we show the patterns obtained from ray 
dynamics for both M-type [Fig. 4(a-b)] and I-type [Fig. 
4(c-d)] functions for different values of n. 
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FIG. 4 Ray-optics prediction p(x, y) of lasing emission for 
the M-type (a-b) and I-type (c-d) wave functions, where (a) 
n = 3.33, (b) n = 1.85, (c) n = 6, and (d) n = 3.7. The thick 
black lines in (b-d) means silvered surfaces. 

Consider the M-type QBS [Fig. 2(a)] that lives on the 
small islands of the PMs of Fig. 1, as shown in [19]. Let 
us also assume that the flat surface of the cavity is sil- 
vered, so light bouncing inside the cavity may only escape 
through its deformed boundary. In this case the M-type 
QBS is equivalent to the bow-tie mode present in the 
ARCs with large deformations [16]. As an illustration, if 
the waveguide is constructed with a semiconductor with 
refraction index n = 3.33, the line for refractive escape, 
shown in Fig. 1(a), will cut the small islands centered 



at sin(x) = 0.3 approximately through the middle. Since 
the average angle of incidence of the M-type function on 
the deformed boundary is slightly less than the critical 
angle Xc ~ 17.5°, the reflectivity R of this boundary is 
(given by the Fresnel's coefficients) sufficiently high (~ 70 
percent) to permit laser action with the direction shown 
in Fig. 4(a). In contrast, if light is allowed to escape only 
through the flat boundary then the M-type function has 
two different angles of incidence on that boundary, see 
Fig. 1(b): one at Xc ~ 90° (the two small half- islands 
at sin(x) w 0) and other at Xc ~ 33° (the small island 
centered at sin(x) w 0.57). For a semiconductor cavity 
with n < 6, the boundary presents a low R (< 50 per- 
cent) at normal incidence (x c ~ 90°). To reduce the laser 
threshold we could leave just an unsilvered region arround 
x = L/2 to permit the escape of light only for Xc ~ 33°. 
In this situation a refraction index of n = 1.85 leads to 
a reflectivity of ~ 70 percent. The emission under these 
conditions is shown in Fig. 4(b). 

The I-type QBS of Fig. 2(b), living in the large is- 
land centered at x — L/2 and with a small sin(x) in both 
PMs of Fig. 1 (as shown in [21]), can be directly asso- 
ciated with a curved mirror Fabry-Perot resonator since 
it bounces with almost normal incidence between both 
walls of the cavity. We will consider only the flat bound- 
ary of the cavity to be completely silvered, but the next 
arguments also work in the opposite case. To construct 
a kind of Fabry-Perot resonator with "one semi-silvered 
mirror" we need a semiconductor material with n ~ 6 in 
order to have a reflectivity of about 50 percent (R ~ 0.5), 
see Fig. 4(c). But note that since n is quite large, the 
beam spreads out due to Snell's law. This problem may 
be fixed by also silvering the deformed boundary of the 
cavity but leaving a hole smaller than the width (in real 
space) of the I-type resonance. For example, if the unsil- 
vered hole is half of the width of the I-type function, the 
reflectivity will increase up to 65 percent. However, under 
this condition, already a value of n = 3.7 produces a 50 
percent reflectivity which corresponds to a Fabry-Perot 
resonator with a semi-silvered mirror, see Fig. 4(d). 

In Rcf. [15,16] it was shown that the ray dynamics 
prediction of lasing emission is in rather good agreement 
with experimental realizations, so we expect that the 
predictions of Fig. 4 are also valid. 

It is clear that by an appropiate choice of the geomet- 
rical parameters of the cavity, one can tailor design the 
emission pattern. For example, V-type QBS may be pro- 
duced with (d,a,L) — (1.0,0.51,5.55), as shown in Fig. 
5(c). Figs. 5(a) and 5(b) are plots of the PMs for the 
deformed and flat boundaries of the cavity, respectively. 
In Fig. 5(c) the wave function of mode number 18 for 
a resonant energy (supporting 20 propagating modes) is 
shown. This is a V-type QBS living inside the stability 
islands of the PMs of Figs. 5(a-b). As before, if the flat 
surface of the cavity is silvered, the normal incidence of 
the V-type QBS on the deformed boundary will produce a 
kind of Fabry-Perot resonator with a "semi-silvered mir- 
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ror" when using a material with n = 6 (R ~ 0.5), shown 
in Fig. 5(d). Finally, with a deformed silvered boundary, 
the reflectivity will be about 70 percent for an index of 
refraction n = 2.5, see Fig. 5(e). 



and alloys [27]. 
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FIG. 5 Poincare Maps for the closed cavity with surface 
of section (a) at the deformed wall, and (b) at the flat wall. 
(d, a, L) — (1.0,0.51,5.55). Critical lines for refractive scape 
are shown, (c) V-type QBS | ^(x,y) | 2 corresponding to 
left-incomming mode 18, using E = 4095.1021. Ray-optics 
prediction p(x,y) of lasing emission for (d) n = 6 and (e) 
n — 2.5. The thick black line in (e) means a silvered surface. 

The quality factor Q gives a measure of the response 
of the cavity to excitation. Following the procedure 
of Ref. [?], we obtain for our cavity the expresion 
Q = 2TTy/R~nl/(l — R)\, where / is the length of the 
trapped trajectory. The reflectivity R depends on the 
index of refraction n and, as in the case of Fig. 4(d), 
on the size of the unsilvered hole. For I-, V- and M- 
type QBS, I is approximately d + 2a, 2.2(d + 1.57a), and 
4.4(d + 1.45a), respectively. Hence Q is approximately 
four times larger for M-type that for Ftype QBS, for a 
given A <C I. 

In conclusion, we have proposed the experimental real- 
ization of microlasers as open semiconductor resonators 
with controlable laser emission. Using ray as well as wave 
dynamics we explained the physical mechanism by which 
lasing can occur; it takes advantage of mixed chaotic dy- 
namics together with quantum tunneling into classically 
forbidden phase space regions. Finally, it is worth to 
point out that the values of the index of refraction used 
here are typical of conventional semiconductor materials 
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